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Abstract: The double-Kerr solution is generated using both a Backhand transformation 
and the Belinskii-Zakharov inverse-scattering technique. We build a dictionary between o 
the parametrisations naturally obtained in the two methods and show their equivalence. 
We then focus on the asymptotically flat double-Kerr system obeying the axis condition 
which is invariant; for this system there is an exact formula for the force between the 
two black holes, in terms of their physical quantities and the coordinate distance. We then 
show that 1) the angular velocity of the two black holes decreases from the usual Kerr value 
at infinite distance to zero in the touching limit; 2) the extremal limit of the two black holes 
is given by \J\ = cM 2 , where c depends on the distance and varies from one to infinity as 
the distance decreases; 3) for sufficiently large angular momentum the temperature of the 
black holes attains a maximum at a certain finite coordinate distance. All of these results 
are interpreted in terms of the dragging effects of the system. 
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1. Introduction 

Any relativistic gravitational theory which reduces, in an appropriate limit, to Newtonian 
gravity, must contain gravitomagnetic effects [|J. In particular these include forces between 
mass currents and therefore spin-spin forces between classical rotating (and therefore ex- 
tended) objects. In General Relativity, the force on a spinning body at rest in the exterior 
field of an arbitrary, stationary, rotating source is given by, to leading order in the distance 
0, 

-J-J + 3(f • j)(f ■ j) \ J=jfJ=Jr p _ GSjJ^ ^ (L1) 




where j, J are the spins of the test body and the source of the gravitational field, respec- 
tively, and r = rf is the separation vector. This result was obtained by considering a 
spinning test particle || in the linearised Kerr field. It is nicely interpreted as a gravit- 
omagnetic dipole-dipole interaction, with exactly the same form, except for the coupling 
constant and the sign, as its counterpart in magnetostatics. For spins aligned with the 



direction of separation, the result (1.1) yields a repulsive/attractive force for parallel/anti- 
parallel spins, in agreement with general arguments Q. 

A natural question is if a repulsive spin-spin force can balance the attractive grav- 
itational force between two classical spinning sources. A naive estimate, based on two 
spinning sources with spin J, mass M and ( |1 . 1| ) reveals only that, if this is to be achieved, 
their separation will be of the order of their minimum size, J/Mc. The question cannot 



be answered at the level of a test particle approximation, since such balance implies that 
the approximation breaks down 0] . To avoid the complexity of an exact solution, one can 
consider the interaction between two spinning particles by an approximation method: a so- 
lution to Einstein's equations up to quadratic order in the metric perturbation around flat 
space, i.e. a post-post Newtonian approximation. Imposing asymptotic flatness one finds 
the suggestive result for the force between the two spinning particles (with spins aligned 
with the direction of separation) [Q] 

Fa _ a MM + a^ i (12) 

where Mi, Ji are the masses and angular momenta of the two particles and £ is the coor- 
dinate distance (in Weyl canonical coordinates) between them. This force, computed from 
the energy momentum tensor associated to the conical singularity at the axis between the 
two particles, is just the sum of the Newtonian and the spin-spin term that one might have 
naively expected, confirming the form ( |1.1|) for the latter. However, ( |1.2|) is only mean- 
ingful if there exists an axis, i.e a set of fixed points of the azimuthal Killing vector field, 
between the two spinning particles, at which the conical singularity is defined. The same 
setup reveals an axis condition, J\jM\ + J2/-M2 = 0, which is incompatible with F = for 
positive mass particles. Thus, already at this level of approximation we can foresee that 
no balance can be achieved. 

One might, however, question the validity of the approximation that leads to (|1.2|). 
Firstly, because it is unphysical to consider pointlike classical spinning particles. Secondly, 
because one might suspect that more terms, other than the quadratic ones in M, J given in 
( |1.2| ), might contribute at the same distance approximation. If one wishes to describe the 
interactions between two classical spinning objects by a regular (on and outside an event 
horizon) exact solution of General Relativity, one is led to consider a solution describing 
two Kerr black holes in equilibrium. The double-Kerr solution || was originally obtained 
by a solution generating technique which explores internal symmetries of the Ernst equa- 
tion H governing stationary, axi-symmetric solutions of the vacuum Einstein equations - 
the Backhand transformation Subsequent analysis (see Q and references therein) has 
revealed that, for two under-extreme objects, the solution cannot be made free of singu- 
larities and hence force balance cannot be achieved. The force between the two Kerr black 
holes takes the form, in a large distance expansion, fl 1 



G 2 (M x M 2 y 3 / Ji h 
c 4 U C J c± \M,C + M 2 ( 



(1.3) 



where Mi, Ji are the Komar masses and angular momenta of the two black holes and £ 
the coordinate distance. The leading order term in the product of the two Komar angular 
momenta is exactly of the form ( |l.l[) , including the numerical coefficient, and the only 



two quadratic terms in M, J in this expression are exactly the ones in (1.2). However, as 



This formula is obtained imposing the axis condition and asymptotic flatness. We should note that we 
have corrected a sign in the force formula given in 0] (cf. (2.5) below). 



anticipated above, there are more contributions at the same order in the distance as the 
term of the form These include the general relativistic corrections to the Newtonian 

force between two masses, already present in the double-Schwarzschild solution. But there 
are also two other terms depending on the spin. Moreover the overall contribution of the 
three leading spin terms is always repulsive, which is apparently in contradiction with the 
intuition gained from (|1.1|), and minimised if J\jM\ + J2/M2 = 0. It should be noted that 
an exact, meaningful force formula, in terms of Komar masses and angular momenta, for 
the general double-Kerr solution seems extremely hard to obtain. Observe that, for such 
formula to be meaningful, one would need to impose, with full generality, the condition 
of vanishing NUT charge (which is generically present in the solution), besides the axis 
condition. 

The double-Kerr solution can also be obtained using a different solution generating 



technique - the inverse scattering method [10, 11 1. This technique has been greatly explored 



over the last few years in higher dimensional general relativity (see [|T^] for a review). In 
particular, in the method was used to generate a double Myers-Perry black hole, a 
five dimensional analogue of the double Kerr. The latter solution was generated using the 



inverse scattering method in [14]. However, the conclusions of [14] are qualitatively different 
from those of ||, namely that, for the physical case when the NUT charge vanishes, the 
spin-spin interaction is always repulsive but maximised if J\/M\ + J2/-M2 = 0. 

The first purpose of this paper is to clarify the apparent contradiction between the 
results obtained by the two different methods. In section 2 we will derive the Kerr solu- 
tion by both the Backhand and inverse scattering method and build an explicit dictionary 



between the two, showing their consistency. We will conclude that (1.3) holds and correct 



some statements in [14]. 




Figure 1: Representation of the invariant double Kerr. This particular case of the double-Kerr 
solution is asymptotically flat and, for sufficiently low values of the angular momentum J describes 
two equal mass and opposite angular momentum black holes with a strut (conical singularity) at 
the axis between them. This strut deforms the geometry of the horizon, deformation which we have 
made no attempt to represent in the figure. 



The second purpose of this paper is to consider the special case when the double Kerr 
system is invariant (cf. section 3) - figure [l|. The solution can then be reduced to a 
three parameter family (of two asymptotically flat Kerr black holes with an axis between 



them) from which some simple exact results can be obtained. It was already shown in [15 
that the exact force between the black holes is given by 



M 2 



/ 2MG\ 



-i 

(1-4) 



where both black holes have Komar mass M but opposite Komar angular momenta J and 
—J; C is again the coordinate distance (in Weyl canonical coordinates). In this special 
case, the explicit angular momentum dependence vanishes from the force formula, which 
becomes identical to the one of the static case (see, for instance [fig] ). Herein, we will 
further show that: 

1) the exact angular velocity of the horizon of the two black holes is given by 2 

o 1 = — 3 y 2M n 2 , (i.5) 

which varies from the Kerr value at £ — > +oo to zero at £ — > 2M (corresponding to 
the two horizons touching) - figure |3| (left); 

2) The extremal limit of the black holes (i.e when they have zero temperature) depends 

on the coordinate distance, and is given by 



which gives the usual Kerr ratio at £ — » +oo but tells us that, in the invariant 
double Kerr system, we can have a black hole with arbitrarily large intrinsic angular 
momentum for fixed mass, by decreasing the coordinate distance; 

3) The temperature of the black holes, for sufficiently high angular momentum, attains a 
maximum at a certain finite coordinate distance, cf. figure ^ (right). 

Moreover, we will interpret all of these results, which will be presented in section 3, in 
terms of the mutual dragging of inertial frames between the two black holes. We shall also 
offer a physical interpretation to the three spin interactions in ( |1.3[ ) which is consistent 
with the intuition gained from (1.1). 



We shall make some final remarks in section 4. 



2 Geometrised units, G = 1 = c are used throughout this paper beyond this point. 



2. The Double-Kerr solution 



The Double-Kerr solution is an exact, seven parameter family of vacuum solutions of four 
dimensional General Relativity, describing a stationary, axi-symmetric spacetime contain- 
ing two Kerr black holes (or two naked singularities or one of each) . Physically, the seven 
independent parameters can be thought to be 

M 1 ,M 2 ,J 1 ,J 2 ,(,b NUT ,M axis ; (2.1) 

these are, respectively, the two Komar masses and the two Komar angular momenta of the 
individual black holes, the distance between them, the NUT charge of the spacetime and 
the Komar mass associated to the axis between the two black holes. The latter parameter 
is zero if the axis condition is obeyed. The existence of a Komar mass associated to the 
"axis" has been discussed in jl7|], in the context of the double Kerr solution and in 13], in 
the context of the double Myers-Perry solution. 

We will now derive the solution using two different solution generating techniques. 
The full metric is extremely involved and we shall not present it explicitly, although all the 
ingredients for its construction will be displayed. We will, however, give explicit expressions 
for the relevant physical quantities. 



2.1 Backlund Transformation 

A general stationary, axi-symmetric metric can be written in the form 

ds 2 = -f[dt - oodcf)} 2 + - [p 2 d(f> 2 + e 2 ^{dp 2 + dz 2 )} . 

All metric functions /, u, 7 depend on p, z only. Solutions to the vacuum Einstein equations 
are determined by the Ernst potential £ , which we write in the form suggested by Yamazaki 

m ' 




For the case of interest herein, E± are 2x2 complex matrices defined as 



n j = , cr j = yP 2 + ( z - a j) 2 , 

where i = 1,3 (k = 2, 4) are line (column) indices and an- = a, — a/%. The solution is 
determined by the 3 real independent distances between the four points along the z axis 
a i < a 2 < a 3 < 04 an d the 4 real parameters u\, u)2, (^3, W4. The three metric functions are 
now determined as follows: 

/ = 4 » A = M(E^E* + ) , B = E+E* + ; 



E± = det 



rij + n k 
a-ik 



± 1 



where ' * ' denotes complex conjugation, -ft (Sj) stands for real (imaginary) part and 



L = ^ det 



t=l,3 



(n t - rc fc )d it H (1 - tf it ) 



, M = - ^ det 



t=l,3 



-(a t0 + a fc0 )c)jt H (1 

a-ik 



ajo = ctj — z and is a constant determined by the condition 

lim uf = , 

meaning that 4> = const is non-rotating at infinity (when the NUT charge vanishes); 



P 2 7 _ 

6 ~ k ST> 



S = det 



2^ 



T = det 



2<7fc 



where fc is a constant determined by the condition 

e 2 ^ 

lim — — = 1 , 

p->oo / 

meaning that there are no conical singularities at infinity. 

Let us notice that instead of the seven parameters mentioned above it is more conve- 
nient to work with the following set 



mi,m 2 , Ai, \ 2 ,ai,a 2 ,(, 



(2.2) 



related to the previous set by 



021 ^2 - ^1 

mi = sec 

2 2 



ai 



a 2 



C = 032 + 
^3 + ^4 



021 + «43 



043 L<J4 — 0^3 
m2 = ^" SeC ^ _ 



Ai 



a; 2 - uji 



A 2 



W4 — W3 



2 2 2 2 

We now consider the most relevant quantities to analyse the physics of this system. 
Consider five regions denoted /, II, III, IV, V according to z < a\, a\ < z < 02, a 2 < z < 
03, 03 < z < C14, z > 04 and always p = 0. 

• The condition for asymptotic flatness is that 

lo\i = Ibj^uT = —uj\v = . 

This means that the overall NUT charge, b^UT, vanishes. Explicitly this yields 

a% C 2 + 2mim 2 b-.C + a°_{m\ cos 2 A 2 - m\ cos 2 Ai) = . (2.3) 

We have defined 



a± = mi sinai±m2 sin 02 , b± = cosai sin A2±cosa 2 sin Ai , a ± = n%\ sin Ai±m2 sin A2 . 



The axis condition, which guarantees that p = is a set of fixed points of the 
azimuthal Killing vector field, is 

to\ni = . 

Explicitly 

— [( 2 — m\ — m\+(a x _) 2 +2mim2 cos(a2— an)][o" C 2 +2"i4m2&+C+a+("i| cos2 M—m\ cos 2 Ai)] 

+4mim2 sin(a2 — «i)[C 2 ~~ m i — m \ + ( a -) 2 ]C 
—Am\m2 cos(a2 — ai)[a^_( 2 + a^(m 2 cos 2 A2 — m\ cos 2 Ai)] = . (2.4) 

The force between the two black holes is 

p _ ]_ /j _ e - 7 | /7 \ mim 2 cos (a 2 -«i) ^ ^ 

4 C 2 — m\ — m\ + (a^_) 2 — 2m\m2 cos(a2 — ol\) 

Note that, since this force is computed by the conical singularity in region III, it 
is only meaningful if the axis condition - which guarantees that there is an axis in 
region III, and hence that discussing a conical singularity is meaningful - is obeyed. 
Otherwise this is just a formal expression. 
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It is also relevant to note that the expressions for the Komar masses and angular 
momenta of the individual black holes take particularly simple (algebraic) forms in terms 
of the imaginary part of the Ernst potential \& = Q£ (the notation \P| z=aj assumes implicitly 
that p = 0): 



o 

00 



M 1 = -~u\ II {V\ g=aa -*\ z=ai ) , M 2 = --u\ IV {^\ z=aA -% z=a3 ) , (2.6) 

° 

Ji = \"\n (M 1 - °f) , J 2 = \u\ IV (M 2 - a f) . (2.7) O 

2.2 Inverse Scattering Method 

In four spacetime dimensions, the inverse scattering method (or Belinskii-Zakharov method) 
|1C, |ll]] can be used to construct new Ricci flat metrics with two commuting Killing vector 
fields from known ones, by using purely algebraic manipulations. Such metrics can always 
be written in the form 

ds 2 = G ab (p, z)dx a dx b + e 2u( - p ' z) (dp 2 + dz 2 ) , (2.8) 

where a, b = 1, 2. To use the BZ method, we need a seed solution, which we take to be the 
double-Schwarzschild solution |l9|, |(]] . It is defined by 

G = diag<^ , p } , (2.9) 

I P2P4 MiA*3 J 

e 2u = ^P2^(piP2 + p 2 ) 2 {p2V3 + p 2 ) 2 {piP4 + p 2 ) 2 {p3Pi + p 2 ) 2 ^ ^ 

1*11*3(1*11*3 + p 2 ?{^ + p 2 ) 2 (iA + p 2 ){i4 + p 2 ){i4 + p 2 )(pI + p 2 ) ' 



-7- 



(-(g+A/),l) (-(g+A m ),l) (-(g+Av),l) 


Figure 2: Rod structure for the double Kerr spacetime. Next to each rod the corresponding 
eigenvector is displayed. For the double Schwarzschild solution these are simply (1,0) for the 
timelike rods and (0, 1) for the spacelike ones. 

where A: is a constant to be determined below and 

Pk = vV + ( z ~ a k) 2 ~ (z - a k ) , p k = -\J p 2 + (z - a k ) 2 - (z - a k ) . 
Its rod structure is represented in figure 



Using the method suggested by Pomeransky [22], we remove two anti-solitons, at z = 
01,03 and two solitons, at z = 02,04. Thus, the seed metric becomes conformal to flat 
space 

G = diag{-l,/^} = G . (2.11) 

MlA*3 Ml^3 

We will actually take the rescaled metric Go to be our seed (bearing in mind that one 
should multiply the final metric by the overall factor P2P4/ HiPs)- We take the generating 
matrix to be 

§o(A, p, z) = diag {-1, p 2 - 2zX - A 2 } . (2.12) 

One can verify that this matrix solves the Lax pair constructed in the BZ method (see 
|1C. |ll]]). The double Kerr black hole is now obtained by a 4-soliton transformation: using 
Go as seed, we add: 

• one anti-soliton, at z = a\ with BZ vector = (1,20464); 

• one anti-soliton, at z = 03 with BZ vector m Qb ' = (1,20363); 

(2) 

• one soliton at z = 02 with BZ vector m Qb = (1, 20202); 

• one soliton at z = 04 with BZ vector = (1, 20404). 

Notice that we have introduced four new parameters: 61, 63, C2, C4. 

Following the BZ algorithm one finds the double-Kerr solution. The t, (ft metric is given 

by 

to** » A (Go) ac 4 fc) (f-) H 4°(GoU 



— ; <^ab — {<^0)ab ~ / J 



^ 1 PkPi 



where p k = p k for k = 2,4 whereas p k = p k for k = 1,3. The four vectors nv® have 
spacetime components: 

^ l (fi k/ ,p,z)} , (2.13) 



and the four BZ vectors where given above. The symmetric matrix T, whose inverse is 
r _1 , is given by 

f rn { a\Go) ab m { P 

i-kl = 2 , - • K Z - l V 

The final quantity we need to have the complete metric is v, which is given by 

e 2u = e 2u detTki 

where r^ ) and T are constructed as in ( |2.14| ) using Go and G, respectively. 

The rod structure of the double-Kerr solution, depicted in figure is completely 
defined by the quantities 



A, 



Gt<f> 



where ao = — oo, 05 = +00. Note that the metric G is the one obtained directly from the 
BZ algorithm. Explicitly 

^ _ 2 (QglQjjjl + Q32Q43^3)Q42 + («41«43C2 + «2lQ32C4)Q31^1^3 

031042(1 + O1C263C4) + 032041(6102 + o 3 c 4 ) + a 2 ia 43 (o 3 c 2 + b 1 c 4 ) 
(021041 - o 3 2a 43 6i6 3 )a 4 2 + (a 4 ia 43 c 2 + a 2 ia 32 c 4 )a 3 i6 3 



A// = 2 



A/// = 2 



Ajv = 2 



Ay = 2- 



031042(01 - o 3 c 2 c 4 ) - a 32 a 4 i(c 2 - &i 6 3 c 4 ) + a 2 ia 43 (6i6 3 C2 - c 4 ) 

(a 2 lQ42C2 ~ Q3lQ4303)Q41 + (02lQ3lC4 ~ 042043^1 )Q3203C2 

032041(1 + 61C263C4) - a 2 ia 43 (6i6 3 + c 2 c 4 ) + a 3 ia42(oic 2 + O3C4) 

(Q3IQ4I + Q32Q42PlC2)043 + (o4lQ4203 ~ Q3lQ32C4)Q2lC2 

021043(01 - 6 3 c 2 c 4 ) + a 32 a 4 i(6 3 - 0ic 2 c 4 ) + a 3 ia 42 (6i6 3 c 2 - c 4 ) 
(a 2 ia 32 C2 + a 4 ia 43 c 4 )a 3 i + (a 32 a 43 &i + a2ia4i03)o 42 c 2 c 4 



03ia 42 (l + 6ic 2 o 3 c 4 ) + a 32 a 4 i(6ic 2 + o 3 c 4 ) + a 2 ia 43 (o3C 2 + 61 c 4 ) 
Moreover, g turns out to be simply expressed as 3 

Xi + Xv 
q = — ' 

These quantities define: 

i) the angular velocity of the horizon of the first and second black holes: 



q + A// "- q + Ajy 



3 q could be defined by the coordinate transformation dt — > dt — qd(f> which ensures that the resulting 
metric obeys linip^oo Gt<j> = 0, when the NUT charge vanishes. 



ii) the NUT charge of the spacetime b^uT] in particular, asymptotic flatness holds iff 



-(q + X I ) = 2b NUT = q + X v = , (2.16) 

iii) the axis condition 

q + A m = , (2.17) 

iv) the regularity condition; choosing a periodicity for (j) such that the only potential 

conical singularity is at p = 0, 02 < z < (13, the regularity condition turns out to be 

denominator[q] = denominator [Xin] ■ (2-18) 

2.3 Dictionary 

There are seven independent parameters in the double Kerr solution. The ones that natu- 
rally arise in the inverse scattering method are 

021, a 32 , «43,&i ; h, C2, c 4 • (2.19) 
These are related to the set that arose in the Backlund transformation (EO) via 



a 21 v / (l + 6f)(l + cj) a 43 V(l + ggj + 4) 

flJi = - — , 77l 2 = - -, — 

2 l + 6ic 2 2 1 + 6 3 C4 

, \ c 2 - b\ c 4 - 63 

tanAi = : , tanA 2 



1 + 6lC 2 ' 1 + 63C4 

c 2 + bi c 4 + 63 

tanai = : , tana2 



1 - 61 c 2 ' 1 - 63C4 

a 2 i + a 43 
C = °32 H 2 ■ 

Inverting these relations we have 

<x 2 i = 2mi cos Ai , 032 = C — m i cos -^l — m 2 cos A2 , a 4 3 = 2m 2 cos A 2 , 

ai - Ai ai + Ai a 2 - A 2 a 2 + A 2 
01 = tan , c 2 = tan , 03 = tan , c 4 = tan . 

It is now a matter of algebra to verify that, using this dictionary, ( |2.3[ ) and ( |2.4[ ) exactly 
match ( p. 16 ) and ( 2.17] ). Moreover (2Tq) is equivalent to F = in (|2.5|), i.e. 



mim 2 cos(a 2 — a\) = 



The latter equation together with ( |2.3| ) and ( |2.4| ) can be verified to exactly match expres- 
sions (4.176), (4.186) and (4.20) in H. Thus the results obtained with the inverse scattering 
method and Backlund transformation are consistent, as they should. 



4 Choosing the periodicity of in this way amounts to choosing the constant k such that linip^oo e 2v — 1. 
The calculation of the regularity condition amounts to computing the period of <f> in region III, which can 
be done using the procedure described in j2l]]. Again note that if the axis condition is not obeyed, the 
result is formal, with no meaningful interpretation. 



We can now comment on the results of [14] mentioned in the introduction. The dis- 
crepancy originates from the fact that the quantities identified therein as being the black 
holes mass and angular momentum are not the Komar quantities of the black holes. Thus, 
although the force formula written in [jl4j (formula (34) therein) is correct, and indeed 
equivalent to formula (2.5) above, its interpretation is not. 5 Moreover, the force formula 
is only meaningful, as we have emphasised, when both the axis condition and asymptotic 
flatness are guaranteed. But the condition imposed in 14] to guarantee vanishing NUT 
charge (t>i = 0) does not guarantee asymptotic flatness, as follows straightforwardly from 
(|2~20|) and Q . 



2.4 Physical Quantities 

The calculation of the Komar masses and angular momenta of the individual black holes 
can actually be simplified by using both methods. Noting that 



uii = -(q + A//) 
then (|2~6| ) and ([T]]) become 
7 + Xn 



uiv = -(<? + A/y) , 



Mi 



z=a 2 



z=a\j i 



M 2 = i±^K(* 



Z=C14 



z=a 3 



+ A// 



M!-^l 

2 



Jo 



q + A/y 



(2.21) 



(2.22) 



The point is that whereas the ^'s are simpler to compute using the Backhand method, the 
(g + A)'s are simpler to compute using the inverse scattering method. The former are 



«32C2 + Q43C4 + Q42&3C2C4 

a 42 + a 43 6 3 c 2 + a 32 b 3 c 4 

021^1 + 04204 + CI4161C2C4 



z=a2 



a 3 ibi + a 43 C4 + 041616304 

(J4I + 043^1^3 + A31&3C4 

a 2 i6 2 + 03203 + 031^1 ° 3 c 2 



(J41 + a42^ic 2 + a 2 ic 2 C4 ' " " : a 3 i + a 32 6ic 2 + a 2 i6 3 c 2 

We shall now observe a couple of consistency checks on the physical quantities com- 



puted. Firstly, for each black hole a Smarr formula should hold 

Ai Tj 



Mi 



+ 20f Ji 



(2.23) 



The of each black hole can be shown to be 

031042(01 - o 3 c 2 c 4 ) - a 32 a 4 i(c 2 - hb 3 c 4 ) - a 2 ia 43 (c 4 - hb 3 c 2 ) 



A x = 2ira 2 i(q + X H ) 
A 2 = 27ra 43 (g + A/y) 



o 3 ia 42 (l + M3C2C4) + 032041 (bic 2 + 6 3 c 4 ) + a 2 ia 43 (6ic 4 + b 3 c 2 ) 

«21«43(frl ~ O3C2C4) + «32Q4l(03 ~ OlC 2 C 4 ) - a 3 ia 4 2(c4 ~ Ol^) 

031042(1 + &1O3C2C4) + a 32 a 4 i(6iC2 + 6 3 c 4 ) + a 2 ia 43 (6ic 4 + b 3 c 2 ) 



5 We can easily build a dictionary between the parametrisation used in [jwj, which we denote a,i,bi,fhi,Zi, 



and the parametrisation used in section 



2.1 



rrii = rrii cos a. 



bi — rrii sin a. 



a* = mi sin Ai 



Z2 — Zi — £ 



(2.20) 



Moreover, using the method described in [21] to compute the temperatures, these can be 
shown to be very simply related to the areas: 



021 



_ T2= M- 



It follows, using also the results in section 2.2 that the Smarr formula ( 2.23| ) is obeyed for 
each black hole. 

Secondly, we would like to check that, when b^UT = 0, the ADM mass/angular mo- 
mentum is the sum of the Komar masses/angular momenta of the individual constituents 
of the spacetime. Note that we say "constituents" rather than black holes because, as 
pointed out in [p^[ , when the axis condition is not obeyed, the rod in region III has a 
non-vanishing Komar mass and angular momentum. These are given by 



+ Am 



z=a3 



^1 



l) 



J n.T.i 



q + A/// 



M nxi 



Q32 
2 



(2.24) 



\z=a2) 1 

It can thus be verified that, for b^uT = 0, 
M ADM = M 1 + M axis + M 2 

_ 031042(021 + 043)(1 ~ M3C2C4) + 032041 («43 ~ Q2l)(ftlC2 ~ b 3 C A ) + Q2ia43(a41 + a32){b3C2 

2(031042(1 + 6163C2C4) + a 32 a 4 i(6iC2 + 63C4) + 021043(6302 + 61C4)) 

and a similar expression, albeit more involved, holds for J ad Ad- 
Let us conclude this section explaining how to obtain formula (1.3) 0. One introduces 
two dimensionless parameters e-i = rrii/C,. Then one imposes the axis and asymptotic 
flatness conditions to determine ai, 02 in terms of the other parameters to quadratic order 
in ej. These expressions for the a^'s are then implemented in the formulae for the force, the 
Komar masses and angular momenta. The resulting expressions (to quadratic order) can 
be used to write the force in terms of the Komar quantities and the coordinate distance. 
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3. The Double-Kerr solution with 7^ symmetry 

To interpret the exact force formula (|2.5| ) it should be written in terms of the physical 
quantities ( |2.1[) . Moreover, a meaningful interpretation would require that the conditions 

M axis = = bNUT , (3.1) 

are obeyed. Although it does not seem possible to solve these conditions with full generality 
0, there is a special case of particular interest: when the system has a iJ^ invariance. By 
this we mean a variation of the usual Z2 invariance (z — > —z) that also reverts the <f> 
direction: in the coordinates fl2.8|) the Z2 is a discrete transformation that acts as 

z-*-z, 0-^-0. (3.2) 

Physically, invariance under ( |3.2| ), means that the system is invariant under a it rotation 
along an axis orthogonal to the z axis and which intersects the middle point of the strut 
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between the two black holes. To realise this symmetry in the double-Kerr system, we 
require ( |3,1| ) and, in terms of the Komar masses and angular momenta of the black holes 



Mi = M 2 = M , Ji = - J 2 = J . (3.3) 

The resulting system has three independent quantities M, J, £. Its ADM mass and angular 
momentum are Madm = 2M, Jadm = 0. Although the spacetime has no overall angular 
momentum its individual constituents have. 

3.1 Physical Quantities 



In terms of the BZ parameters ( [2.19 ) the constraints ( [3.1| ) and ( [3.3j ) are implemented by 
the intuitive conditions 

G&21 = «43 = a , h = C4 = b , C2 = 63 = c . 

The axis condition is immediately obeyed with no further restrictions. Asymptotic flatness 
determines function of b, c, £ 

(c + fe)(l + fec) 

a ~ (c-b)(i-b c y ' 

The Komar mass and angular momentum can also be written in terms of b, c, ( as 

~ \c-b) 2 ' J ~ 2(1 -6c) \c-b) C • 

Note that the values of b, c are restricted by the condition that £ > a > and M > 0. The 
trivial case (no black holes) is b = — c. The static case (double-Schwarzschild) is obtained 
by taking 

fe,c^0, -e]-l,0[; (3.4) 

c 

it degenerates into flat space and a single Schwarzschild in the two limits of this interval, 
respectively. The extremal limit (zero temperature) is 1 + be = 0, which follows from the 
temperatures 

1 ( C - 6)2(1 + be) 
1 2vr 4c(c + 6)C 2 ' 

Note that the black holes have the same temperature. However they are not in thermal 
equilibrium since they have a different chemical potential (symmetric angular velocities). 6 
Note also that the extremal limit does not correspond to J = M 2 , as for a single Kerr 
black hole, but rather J = cM 2 . Taking further c = 1 (i.e. independent of C) one ends up 
with the flat space limit. The angular velocity of the horizons are 



n ?=-HSii~ nf 



°We would like to thank R. Emparan for this observation. A similar situation was discussed in 



It can be verified that, in the allowed regions of the b, c parameter space, Jj has the same 
sign as Op, as expected. The limit b — > (c / 0) corresponds to a = £, i.e. the two 
horizons touch. In this limit the force 

F _ (c + b) 2 
166c ' 

diverges and the angular velocity of the horizons vanishes. We shall interpret the latter 
conclusion in the next section. 

Re-expressing the last three formulas in terms of M, J, £, we obtain 



M 2 



2M\ 2 



which is (|1.4p and 



n 1 = 2 -^ = -n 2 , Tl = ± M f M - a \=T 2 , (3.5) 

4J ' 8tt2J 2 (C-2M) ' V ' 



where 



/ „ J 2 C-2M . . 

a = 2\M 2 -—.\ — , 3.6 

y M 2 ( + 2M y ' 

which in particular yields (|1.5|). Note that the extremal limit is a — > 0, which yields (|1.6j), 
and only corresponds to J 2 = M 4 for Q — > oo. The area of the black holes is 

C M(2M-a) V ; 
Note that, in the large distance limit ^ — > oo, 



2M(M 2 + VM 4 - J 2 ) ~ 4vrM (m 2 + x/M 3 - J 2 



A x = 8vr [M 2 + \/M 4 - J 2 ) = A 



which are the standard quantities for a single Kerr black hole. One also observes that in 
the touching limit £ — > 2M all quantities become independent of the angular momentum: 

Oi,O 2 ^0, Ti 



16vrM 
Ai = 32vrM 2 = A 2 . 

Let us remark that the touching limit becomes a coalescence limit corresponding to a single 



Schwarzschild black hole if one takes simultaneously the limit (3^). This follows from the 
fact that the rod structure displayed in figure |2| becomes, in such case, the one of a single 
Schwarzschild black hole. 




3.2 Physical interpretation 

Let us start by analysing the angular velocity, temperature and area for the l!^ invariant 
double Kerr. These quantities are displayed in figures [3| and |j. 

From figure || (left) we see that, as we decrease the coordinate distance, the angular 
velocity of the two Kerr black holes decreases monotonically from the Kerr value (at Q — > oo) 
to zero in the touching limit (C —* 2M). This can be attributed to the "dragging of inertial 
frames" from the other black hole. Take, say, black hole 1. It has a positive intrinsic 
angular momentum and hence a positive angular velocity. But it is dragged in the negative 
direction by the negative intrinsic angular momentum of black hole 2 and therefore it slows 
down, i.e. its angular velocity is smaller than that of an isolated black hole with the same 
mass and angular momentum. And (|ll]) shows that the dragging effect tends to cancel the 
positive angular velocity due to the intrinsic angular momentum as the black holes touch. 
This could have been anticipated on physical grounds; a non-vanishing angular velocity in 
the touching limit would create a discontinuity in the spacetime dragging effects. Let us 



note that it had already been pointed out in [15| that in the double Kerr system the 
black holes would slow down each other, but without relating this fact to the dragging of 
inertial frames of General Relativity. 

The variation of the temperature with the coordinate distance is, for sufficiently low 
angular momentum, essentially the same as in the static case (figure ^ (right), J = 0.2, 0.5). 
In the static case, as we decrease C the black holes become larger (i.e the area increases 
monotonically, which also holds for the stationary case, cf. figure |] (left)) and hence the 
temperature decreases monotonically. However, for sufficiently large angular momentum 
we see a new phenomenon occurring: there is a maximum temperature at a certain, finite 
coordinate distance. We interpret the existence of such maximum as the result of two 



distinct effects. The first one is the one just discussed which also holds for the static case. 
The second one is intrinsically associated to rotation and to the fact that the extremality 
condition, a = 0, depends on the coordinate distance, cf. Q3.6Q . 

Recall that for a Kerr black hole, the temperature decreases monotonically from the 
Schwarzschild limit to the extremal limit. Hence, the temperature can be faced as a mea- 
sure of the deviation from extremality. On the other hand, extremal rotating (uncharged) 
black holes are systems under the maximum allowed mechanical stress. A higher angular 
momentum (and hence higher angular velocity) would destabilise the horizon, transform- 
ing the system into a naked singularity. In the Z2 double Kerr system with fixed M, J 
decreasing the coordinate distance slows down the black holes, as discussed above, due to 
their mutual dragging effects. Thus, if the black holes were initially extremal (and hence 
in their mechanical limit) they relax and become non-extremal as they are brought closer 
- figure U (right) J = 1. Since the temperature is a measure of the deviation from ex- 
tremality, the temperature of the black holes should increase as they are brought together, 
which is exactly what we can see in figure ^ (right), J = 1, for sufficiently large distance. 
As the black holes are brought sufficiently close, they slow down sufficiently so that the 
dominating effect is the same as in the static case, i.e. the increasing in size of the hole, 
and hence the temperature starts decreasing. 




Figure 4: Left: A\ = A2 for M = 1 and various values of J, as a function of the coordinate 
distance £. Right: Parameter a (3.6) measuring deviation from extremality, as a function of £. 



The fact that the extremality limit depends on £ originates a novel phenomenon in 
four dimensional black hole physics. A Kerr black hole can have arbitrarily large angular 
momentum for fixed mass in the Z2 invariant, asymptotically flat (and obeying the axis 
condition) double-Kerr system. This can be achieved simply by decreasing sufficiently £ for 
the given value of the mass. Note, however, that these are not ultra- spinning Kerr black 
holes, since they do not have an arbitrarily large angular velocity. It is well known that 
Myers-Perry black holes (with rotation in a single plane) [24[ have no upper bound for the 



angular momentum per unit mass in spacetime dimension D > 6. In the ultra-spinning 



limit they become "pancake- like" and it has been argued they should be unstable [25]. In 
the present case, by contrast, since the angular velocity is always bounded and actually 
decreases with decreasing £, for fixed mass and for the extremal solutions, we expect the 
horizon not to be too flattened by the rotation. It would, indeed be quite interesting to 
study in detail the horizon (and ergosphere) geometry of these solutions. 

Let us now reconsider the generic double Kerr system, the force formula Ql.3] ) and, 
in particular, the three spin dependent terms. The first observation is that one of the 
Jf, Jf terms will be present if a single black hole has intrinsic angular momentum. Thus, 
unlike the J\ J2 term it is not an interaction between the intrinsic rotation of both black 
holes. Secondly one notes that these J\ , J| terms are always repulsive. These two facts 
suggest that these terms are due to the interaction between the intrinsic spin of one black 
hole and the angular velocity that such intrinsic spin induces in the other black hole, by 
dragging of inertial frames. Thus, the effect is present even if only one of the black holes 
has intrinsic spin and is repulsive, for the dragging originates an induced angular velocity 
parallel to the original spin. If this interpretation is correct one concludes that the spin- 
spin interaction being always repulsive means that the interaction between induced and 
intrinsic spins always puts an upper bound to the interaction between intrinsic spins. 8 

4. Final Remarks 

In this paper we have derived the double Kerr solution using two different solution generat- 
ing techniques: the Backhand transformation and the inverse scattering method. We built 
a dictionary between the parametrisations that naturally arise in either method and used 
it to verify that some of the key quantities to analyse the solution - namely the asymp- 
totic flatness condition, the axis condition and the force between the black holes - match. 
This allowed us to clarify an apparent contradiction between the conclusions of [H] and [14] 



concerning the spin-spin interaction between the two Kerr black holes. Both results are 
technically correct but the interpretation of the result given in the latter paper is not, as 
described in section 2. 

The analysis of the general solution led us to recognise one special case, which we 
have called the invariant double-Kerr, which simplifies sufficiently so that some exact 
statements can be made in terms of physically meaningful quantities, as discussed in section 
3. Of particular interest are the physical consequences of the dragging of inertial frames. 
This is normally discussed in the context of a test particle approximation. It is quite 
striking to clearly see this dragging at the level of an exact solution. One should mention 
this has primarily been achieved in the black Saturn solution p^| , which has the advantage 
of being everywhere regular (on and outside the horizon), but the disadvantage of being a 
five dimensional solution and therefore perhaps less physical. 

7 The fact that the extremality limit has \ J\ > M 2 has a counterpart for the charged dihole ]2(j, |27], Q 
wherein the extremal limit has \Q\ > M. We thank R. Emparan for this observation. 

8 Note that what is meant by 'the interaction being always repulsive' is that turning on the angular 
momentum, for fixed coordinate distance, the force between the black holes always decreases. The same 
does not apply for fixed physical distance. 



Let us conclude with two questions. Firstly, is there a way to understand the area in- 
crease in the invariant system, for fixed mass and angular momentum, as the coordinate 
distance decreases? Note that this holds even for J = and is certainly a consequence 
of the horizon deformation due to the strut. Secondly, one could imagine the following 
"gedankenexperiment" : set two Kerr black holes, with the same mass M and counter- 
aligned spins, J and —J, some distance apart (. Let 



thus, if the conclusions from the double Kerr system are applicable, the two Kerr black 
holes are under-extreme at their position but would be over-extreme at infinite distance. 
Could we give them initial velocities — v and v along the z direction such that they reach 
a distance where they become over- extreme and hence naked singularities? That would 
seem to violate both the second law of thermodynamics and cosmic censorship. 

Note Added in the Revised Version 

After the pre-print of this paper appeared, another paper (arXiv:0809.2422 [gr-qc]) ap- 
peared considering the Zg system. It correctly pointed out that the explicit metric written 
in our previous version (eq. (3.8) therein) did not describe the extremal limit of the Zg 
system, asymptotically flat and obeying the axis condition. The corrected version of such 
metric is, however, much more involved, and therefore we chose not to include it in this re- 
vised version, but rather analyse it elsewhere. Note that this does not alter the correctness 
of any of the other results in this paper. Note also that the original physical statement that 
in the 7L\ system one can have under-extreme black holes with J > M 2 was also stated in 
arXiv:0809.2422 [gr-qc] but, unfortunately, not acknowledged therein that it had been first 
discussed in our paper. 
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